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Abstract. Models developed within the framework of the relativistic impulse 
approximation for quasielastic (QE) electron and neutrino-nucleus scattering 
are discussed. Different descriptions of final-state interactions (FSI) in the in- 
clusive scattering are compared: the relativistic Green's function (RGF) and 
the relativistic mean field (RMF). The results of the models are compared with 
the recently measured double-differential charged-current QE (CCQE) neutrino 
MiniBooNE cross sections. 



1 Introduction 

Several decades of experimental and theoretical work on electron scattering 
have provided a wealth of detailed information on nuclear structure and dynam- 
ics |Hl|2l. Additional information is available from neutrino-nucleus scattering. 
Recently, muon charged-current (CC) neutrino-nucleus double differential cross 
sections have been measured by the MiniBooNE collaboration J3). High-quality 
predictions for neutrino-nucleus cross sections are needed for use in ongoing ex- 
perimental studies of neutrino oscillations at GeV energies and a proper analysis 
of data requires that the nuclear response to neutrino interactions is well under 
control and that the unavoidable uncertainties on nuclear effects are reduced as 
much as possible. Although the two situations are different, electron scattering 
is the best available guide to determine the predictive power of a nuclear model. 
Nonrelativistic and relativistic models have been used to describe nuclear effects 
with different approximations. Relativity is however important at all energies, 
in particular at high energies, and in the energy regime of many neutrino exper- 
iments a fully relativistic approach is required, where not only relativistic kine- 
matics is considered, but also nuclear dynamics and current operators should be 
described within a relativistic framework. 

Models for the QE exclusive and inclusive electron and neutrino scattering 
are presented in this contribution. In the QE region the nuclear response is dom- 
inated by the mechanism of one-nucleon knockout, where the probe interacts 
with a quasifree nucleon that is emitted from the nucleus with a direct one-step 
mechanism and the remaining nucleons are spectators. In the exclusive (e, e'p) 
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reaction the outgoing nucleon is detected in coincidence with the scattered elec- 
tron, the residual nucleus is left in a specific discrete eigenstate and the final state 
is completely specified. In the inclusive (e, e') scattering only the scattered elec- 
tron is detected, the final nuclear state is not determined, and the cross section 
includes all the available final nuclear states. The inclusive CC scattering where 
only the charged lepton is detected can be treated with the same models used for 
the inclusive (e, e') reaction. 

For all these processes the cross section is obtained in the one-boson ex- 
change approximation from the contraction between the lepton tensor, that under 
the assumption of the plane-wave approximation for the initial and the final lep- 
ton wave functions depends only on the lepton kinematics, and the hadron tensor 
W^ v , whose components are given by bilinear products of the matrix elements 
of the nuclear current J M between the initial and final nuclear states, i.e., 

W'» = Y t (*f I J ^(l) I *»•)<*< I J U H<l) I 9f)5(Ei+u-E,), (1) 
/ 

where to and q are the energy and momentum transfer, respectively. Different 
but consistent models to calculate the components of the hadron tensor in QE 
electron and neutrino-nucleus scattering are outlined in the next sections. 

2 The exclusive (e, e'p) reaction 

For the exclusive (e, e'p) reaction a model based on the distorted-wave impulse 
approximation (DWIA) has been developed to calculate the matrix elements in 
Eq. ([T). The model is based on the following assumptions fl]|2l[4]: 

i) An exclusive process is considered, where the residual nucleus is left in a dis- 
crete eigenstate n of its Hamiltonian. 

ii) The final nuclear state is projected onto the channel subspace spanned by the 
vectors corresponding to a nucleon, at a given position, and the residual nucleus 
in the state n. This assumption neglects effects of coupled channels and is justi- 
fied by the considered asymptotic configuration of the final state. 

iii) The (one-body) nuclear-current operator does not connect different channel 
subspaces and also the initial state is projected onto the selected channel sub- 
space. This is the assumption of the direct-knockout mechanism and of the IA. 

The amplitudes of Eq. (fl} are then obtained in a one-body representation as 

K! 2 {X^ I ■?'"(<?) I Vn) j (2) 

where the one-body nuclear current, is the single-particle (s.p.) scat- 
tering state of the emitted nucleon, ip n the overlap between the ground state of 
the target and the final state n, i.e., a s.p. bound state, and the spectroscopic 
factor A„ is the norm of the overlap function, that gives the probability of re- 
moving from the target a nucleon leaving the residual nucleus in the state n. In 
the model the s.p. bound and scattering states are eigenfunctions of a non Her- 
mitian energy dependent Feshbach-type optical potential and of its Hermitian 
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conjugate at different energies. In standard DWIA calculations phenomenologi- 
cal ingredients are usually employed: the scattering states are eigenfunctions of 
a phenomenological optical potential determined through a fit to elastic nucleon- 
nucleus scattering data and the s.p. bound states are obtained from mean-field 
potentials, or can be calculated in a phenomenological Woods-Saxon well. 

The model can be formulated in a similar way within nonrelativistic |4| 
DWIA and relativistic RDWIA frameworks [5]. Both the DWIA and the RD- 
WIA have been quite successful in describing (e, e'p) data in a wide range of 
nuclei and in different kinematics [2,5-0. 

3 Inclusive lepton-nucleus scattering 

In the inclusive scattering where only the outgoing lepton is detected all elas- 
tic and inelastic channels contribute and the complex potential, with imaginary 
terms designed to reproduce just the elastic channel, should be dismissed. Dif- 
ferent approaches have been used to account for FSI. In the approaches based 
on the RDWIA, FSI have been accounted for by purely real potentials. The final 
nucleon state has been evaluated with the real part of the relativistic energy- 
dependent optical potential (rROP), or with the same relativistic mean field po- 
tential considered in describing the initial nucleon state (RMF) [8 ,9 |. However, 
the rROP is unsatisfactory from a theoretical point of view, since it is an energy- 
dependent potential, reflecting the different contribution of open inelastic chan- 
nels for each energy, and under such conditions dispersion relations dictate that 
the potential should have a nonzero imaginary term [ 10 1. On the other hand, the 
RMF model is based on the use of the same potential for bound and scattering 
states. It fulfills the dispersion relations [ 10] and also the continuity equation. 

A different description of FSI makes use of Green's function (GF) tech- 
niques 0J3-[17). In the GF model, under suitable approximations that are basi- 
cally related to the IA, the components of the hadron tensor are written in terms 
of the s.p. optical model Green's function. This result has been derived by argu- 
ments based on multiple scattering theory ifTUl or by means of projection opera- 
tors techniques within nonrelativistic ifTTI and relativistic fT2UT4l frameworks. 
The explicit calculation of the s.p. Green's function can be avoided ifTTlFPil by 
its spectral representation, that is based on a biorthogonal expansion in terms of 
a non Hermitian optical potential H and of its Hermitian conjugate 1-0 . Calcu- 
lations require matrix elements of the same type as the DWIA ones in Eq. (fJJ, 
but involve eigenfunctions of both H and 1-0 , where the imaginary part gives in 
one case absorption and in the other case gain of flux, and in the sum over n the 
total flux is redistributed and conserved. The GF model allows to recover the 
contribution of non-elastic channels starting from the complex optical potential 
that describes elastic nucleon-nucleus scattering data. It provides a consistent 
treatment of FSI in the exclusive and in the inclusive scattering and gives also a 
good description of (e, e') data 111 1111 2l[l8l . 

For the inclusive electron scattering both nonrelativistic GF H 1 1 II 1 51 and rel- 
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ativistic RGF lfl2l models have been considered, for CC neutrino scattering the 
RGF model has been adopted lfl3l . The results of the RMF and RGF models 
have been compared in ifTHl for the inclusive electron scattering and in l20ll for 
the inclusive CC neutrino scattering. An example is shown in Figure [TJ where 
the RGF, RMF, and rROP cross sections of the 12 C(e, e') reaction calculated in 
a kinematics with a fixed value of the incident electron energy (e = 1 GeV), 
and two values of the momentum transfer (q — 500 and 1000 MeV/c) are dis- 
played. Two parameterizations of the ROP have been used for the RGF calcula- 
tions, i.e., the energy-dependent and A-dependent EDAD1 (RGF1) and EDAD2 
(RGF2) fljD. The results of the relativistic plane-wave IA (RPWIA), where FSI 
are neglected, are also shown in the figure. 

The differences between the RMF and RGF results, as well as the differ- 
ences between RGF1 and RGF2, increase with the momentum transfer. At q 
= 500 MeV / c the three results are similar, both in magnitude and shape, larger 
differences are obtained at q = 1000 MeV/c. The shape of the RMF cross sec- 
tion shows an asymmetry, with a tail extending towards higher values of u, that 
is essentially due to the strong energy-independent scalar and vector potentials 
present in the RMF model. The asymmetry towards higher u> is less significant 
but still visible for RGF1 and RGF2, whose cross sections show a similar shape 
but a significant difference in magnitude. At q = 1000 MeV/c both the RGF1 
and RGF2 cross sections are higher than the RMF one in the maximum region, 
but a stronger enhancement is obtained with RGF1 , which at the peak overshoots 
the RMF cross section up to 40% and it is even higher than the RPWIA result. 

The differences between the RGF1 and RGF2 results are consistent with the 
general behavior of the two ROP's and are basically due to their imaginary part. 
The real terms are very similar and the rROP cross sections are not sensitive 
to the parameterization considered. The scalar and vector components of the 
real part of the ROP get smaller with increasing energies and the rROP result 
approaches the RPWIA one for large values of uj. In contrast, the imaginary 
part presents its maximum strength around 500 MeV, being also sensitive to the 
particular ROP parameterization. 

An example for the 12 C , cross section is shown in Figure [2] For 
the RGF model, the RGF1 results are compared with the results obtained with 
the energy-dependent but A-independent EDAI potential (RGF-EDAI). Calcula- 
tions have been carried out with the same incident lepton energy and momentum 
transfer as for the (e, e') cross sections of FigureQ] Also in Figure|2]the shape of 
the RMF cross section shows an asymmetry with a tail extending towards higher 
values of lj (corresponding to lower values of the kinetic energy of the outgoing 
muon TfJ, An asymmetric shape is shown also by the RGF cross sections, while 
no visible asymmetry is given by the RPWIA and rROP results. Also in this case 
the differences between the RGF and rROP cross sections are consistent with the 
general behavior of the ROP's and are due to their different imaginary part. As 
already shown for the (e, e') reaction, the RGF yields in general a larger cross 
section than the RMF. This may reflect the influence of the pionic degrees of 
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freedom, that can be included in a phenomenological way in the imaginary part 
ofmeROPflSEQ). 

The results in Figure [2] present some differences with respect to the corre- 
sponding (e, e') cross sections in FigureQ] In both cases the differences between 
the results of the different models are generally larger for increasing value of the 
momentum transfer. For neutrino scattering, however, such a behavior is less 
evident and clear. In particular, the RGF1 cross section at q = 1000 Mev/c does 
not show the strong enhancement in the region of the maximum shown in Fig- 
ure Q] where the RGF1 result is even larger than the RPWIA one. In the case of 
neutrino scattering the RGF results in the region of the maximum are generally 
larger than the RMF ones, but smaller than the RPWIA cross sections. The nu- 
merical differences between the RGF results for electron and neutrino scattering 
can mainly be ascribed to the combined effects of the weak current, in particular 
its axial term, and the imaginary part of the ROP [ 20 1 . 

4 Scaling functions 

The comparison has been extended to the scaling properties of the different mod- 
els l|18l|20|| . Scaling ideas applied to inclusive QE electron-nucleus scattering 
have been shown to work properly to high accuracy ll22l|23ll . At sufficiently 
high momentum transfer a scaling function is derived dividing the experimen- 
tal (e, e') cross sections by an appropriate single-nucleon cross section. This is 




Figure 1. Differential cross section of the 12 C(e, e') reaction for an incident electron 
energy e = 1 GeV and q = 500 and 1000 MeV/c. Results for RPWIA (dotted), rROP 
(dot-dashed), RGF1 (solid), RGF2 (long dot-dashed), and RMF (dashed) are compared. 
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basically the idea of the IA. If this scaling function depends only upon one kine- 
matical variable, the scaling variable, one has scaling of first kind. If the scaling 
function is roughly the same for all nuclei, one has scaling of second kind. When 
both kinds of scaling are fulfilled, one says that superscaling occurs. An exten- 
sive analysis of electron scattering data has shown that scaling of first kind is 
fulfilled at the left of the QE peak and broken at its right, whereas scaling of sec- 
ond kind is well satisfied at the left of the peak and not so badly violated at its 
right. A phenomenological scaling function f^ xp (ip') nas been extracted from 
data of the longitudinal response in the QE region. The dimensioneless scaling 
variable tp'(q, u>) is extracted from the relativistic Fermi gas (RFG) analysis that 
incorporates the typical momentum scale for the selected nucleus |9 22]. Al- 
though many models based on the IA exhibit superscaling, only a few of them 
are able to reproduce the asymmetric shape of f^ xp (tp') with a significant tail 
extended to high values of oj (large positive values of One of these is the 
RMF model. In contrast, the RPWIA and the rROP lead to symmetrical-shape 
scaling functions which are not in accordance with data analysis [9,24). The 
scaling function of the RGF and RMF are very similar for lower values of the 
momentum transfer (q = 500 — 700 MeV/c) and in good agreement with the 
phenomenological function fl8l . the asymmetric tail of the data and the strength 
at the peak are fairly reproduced by both models, while visible discrepancies 
appear increasing q. 

In Figure[3]the scaling function extracted from the differential cross sections 
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Figure 2. Differential cross section of the 1 C(z/ M , n ) reaction for E„ = 1000 MeV and 
q = 500 MeV/c and 1000 MeV/c. Results for RPWIA (dotted), rROP (dot-dashed) RGF1 
(solid), RGF- EDAI (long dot-dashed), and RMF (dashed) are compared. 
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of the 12 C(v^, reaction evaluated for different models at q = 500 and 1000 
MeV/c (shown in Figure [2]) are compared with the phenomenological scaling 
function extracted from the analysis of (e, e') data. Both the RMF and the RGF 
models give an asymmetric shape, with a tail in the region of positive tp'. In con- 
trast, the RPWIA and the rROP results do not show any significant asymmetry. 
As a general remark, these results for the scaling functions follow similar trends 
to those already applied to the behavior of the cross sections in Figure|2] 

5 Comparison with Charged-Current Quasielastic MiniBooNE data 

The CCQE 12 C(v tll cross sections recently measured by the MiniBooNE 
collaboration Q have raised debate over the role of the theoretical ingredients 
entering the description of the reaction. The experimental cross section is un- 
derestimated by the RFG model, and by other more sophisticated models, unless 
the nucleon axial mass Ma is significantly enlarged (1.35 GeV/c 2 in the RFG) 
with respect to the accepted world average value (1.03 GeV/c 2 1 25 , 26 ]). Before 
drawing conclusions about the value of the axial mass it is however important 
to evaluate carefully the contributions of all the nuclear effects. Within the QE 
kinematic domain the treatment of FSI is essential for the comparison with data. 
The comparison between the results of the RMF and RGF models can be helpful 
for a deeper understanding of the role played by FSI in the analysis of CCQE 
data Il27ll . 

The CCQE double-differential 12 C cross sections averaged over 
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the neutrino flux as a function of 7), for various bins of cos 9, where 9 is the 




Figure 4. Flux-averaged double-differential cross section per target nucleon for the 
CCQE 12 C(i/ M ,^~) reaction calculated in the RMF (blue line), the RGF1 (red), and 
RGF-EDAI (green) and displayed versus T M for various bins of cos 9. In all the calcu- 
lations the standard value of the nucleon axial mass, i.e., M A = 1.03 GeV/c 2 has been 
used. The data are from MiniBooNE (3). 
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Figure 5. Total CCQE cross section per neutron versus the neutrino energy. The cross 
sections calculated in the RMF (blue line), RGF EDAD1 (red), and RGF EDAI (green) 
potentials are compared with the flux unfolded MiniBooNE data of 1 3 1. 
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muon scattering angle, are shown in Figure[4] The RMF results yield reasonable 
agreement with data for small angles and low muon energies, the discrepancy 
becoming larger as 9 and T M increase. The shape followed by the RMF and RGF 
cross sections fits well the slope shown by the data. The two models yield close 
predictions at larger values of T p , for all the bins of cos 9 shown in the figure. 
The RGF cross sections are generally larger than the RMF ones. The differences 
increase approaching the peak region, where the additional strength shown by 
the RGF produces cross sections in reasonable agreement with the data. The 
differences between the RGF-EDAI and RGF-EDAD1 results are enhanced in 
the peak region and are in general of the order of the experimental errors. 

In Figure [5] the total CCQE cross section per neutron obtained in the RMF 
and RGF models are displayed as a function of the neutrino energy and com- 
pared with the "unfolded" experimental data J3|. The RMF results underpredict 
the MiniBoone cross section. Larger cross sections, in particular for larger val- 
ues of E v , are obtained in the RGF with both optical potentials. Visible differ- 
ences are obtained between the RGF-EDAI and the RGF-EDAD1 results, being 
RGF-EDAI in good agreement with the shape and magnitude of the experimen- 
tal cross section and RGF-EDAD1 above RMF but clearly below the data. These 
differences are due to the different imaginary parts of the two ROP's, particularly 
for the energies considered in kinematics with the lowest 9 and the largest T M . 
These kinematics, which were not considered in previous calculations 01311201 . 
give large contributions to the total cross section and emphasize the differences 
between the RGF predictions with the two optical potentials. We notice that 
EDAI is a single-nucleus parameterization, which does have an edge in terms of 
better reproduction of the elastic proton- 12 C phenomenology |fl9l compared to 
ED AD 1 , and also leads to CCQE results in better agreement with data. 

The RMF model generally underpredicts the data [28 1. In contrast, the RGF 
can give cross sections of the same magnitude as the experimental ones without 
the need to increase the standard value of the axial mass. The larger cross sec- 
tions arise in the RGF model from the translation to the inclusive strength of the 
overall effect of inelastic channels. 

These results confirm that before drawing conclusions about the comparison 
with CCQE MiniBoone data and the need to increase the axial mass, the rele- 
vance of all nuclear effects must be investigated. A careful evaluations of the 
relevance of multi-nucleon emission [29, 30] and of some non-nucleonic contri- 
butions ||3T1 would be helpful to clarify the content of the enhancement of the 
CCQE cross sections obtained in the RGF model. A better determination of a 
phenomenological ROP which closely fullfills the dispersion relations deserves 
further investigation. 
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